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Abstract 



In this paper we construct a global, continuous flow of solutions to the Camassa-Holm equation 
on the space H^{M.). In a previous paper 1^, A. Bressan and the author constructed spatially periodic 
f~| ' solutions, whereas in this paper the solutions are defined in all the real line. We introduce a distance 

I functional, defined in terms of an optimal transportation problem, which allows us to study the 

continuous dependance w.r.t. the inital data with a certain decay at infinity. 



1 Introduction 

In [3j the authors present a nonlinear partial differential equation which describes the behaviour of the 
l/^ ' shallow water as a completely integrable hamiltonian system 

^ Ut + 2kUx - Uxxt + ^UUx = "^UxUxx + UUxxx- (1) 

O 



where u is the fluid velocity in the x direction and k is a constant related to the critical shallow-water 
wave speed. For the physical description of such equation, we refer to [HI El EI and to the bibliographic 
references of (2]. In the present paper we study the limit case k = 0, a condition in which, starting 
by a initial smooth data it can develop to a peaked solution with one or more cuspids, the so called 
■ multi-peakon function. As an example, in a simulation shows that starting from a parabolic initial 
^ , data in a periodic domain, the system evolves in a train of positive peakons. The equation can be 
$-H ' written in nonlocal form as a scalar conservation law with an integro-differential source term: 



«*+(y) +Px=0 (2) 



where is defined in term of a convolution: 



1 / 9 U„ 



2 

Observe that the function ^e"'^' is the distributional solution of the equation 

{Id - dxx) f = So 

where Sq is the Dirac measure centered at the origin. A multi-peakon is a function of the form 

N 

I 



=1 



1 



it is well known (see [HllSlini) that if such a function is subject to (HJ, its evolution remains of the same 
shape, and as long as they are well defined, the coefficients pi, qi are the solution to the system of ODE 

N 

i = l,...,iV. 

Pi = Vi'^Pj sign {qi - g-,)e~l'''"^^l, 

In the smooth case, differentiating Q w.r.t. x and multiplying by Ux we obtain the conservation law 
with source term 




The previous equation, together with the Camassa-Holm equation 0, prove that the total energy 

E= [u^{t,x) +ul{t,x)]dx 
Jr 

is a conserved quantity as long as the solution remains regular. Constantin and Escher |H1 Theorem 4.1] 
shown that even if the initial data is sufficiently regular, blow-up of the gradient Ux can occur in finite 
time whenever the initial data has a negative slope. In Section ^ we implement a technique based on 
appropriate rescaled variables in order to resolve the singularities which occurs at the times where the 
gradient Ux blows up. The new system of ODE can be solved in a unique way in a neighborhood of the 
time of blow up and the solution turns out to preserve the energy E also after this time. 

Motivated by the existence of the multi-peakon solutions, whose decay at infinitive is like e"'^', we 
introduce the space Xa of the function with an exponential decay: let < a < 1, we define 

X« = |n G H\R) : = j [u^{x) + ul{x)] e"l^l dx < ooj . (4) 

In this space we define a distance that is related to an optimal transportation problem (see jllj). 
Fetching the theory developed by Bressan and Constantin ^ for the Hunter-Saxton equation, and by 
Bressan and Fonte [2] for the periodic solution of the Camassa-Holm equation, the topology induced by 
functional J constructed in Section El turns out to be weaker than the topology, but useful because 
with this metric we prove the stability of the multi-peakon solutions w.r.t. the initial conditions, as we 
will show in Section [7J 



2 Multipeakon solutions 

In this section we shall construct a solution of the Camassa Holm equation starting from an initial 
condition E Xa of the form 

ul{x) = Y,P^e'\^-'^\. 

The motivation of this choice is given by the form of traveling wave solution (see jH] and Example 
5.2]). Looking for solution of the equation in the traveling wave form u{t,x) = U{x — ct), with 
a function U that vanishes at infinity, the limit of k — > leads to the function U = ce~l^~'^*L The 
evolution of an initial data like Uq remains then of the same shape 

u%t,x) = X^^'iWe"'^'"*^*^' . 
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As long as the classical solution of the problem 



N 



Qi = J2Pj^ '''' '''' 



N 



(5) 



Pi = Pi'^Pj sign {qi -qj)e '''' 



exists, the solution of this system gives the coefficients p{t) = {pi, . . . ,pnJ and q(t) = (gi, . . . , g^r^) for 
the solution u'^[t,x) to the Camassa-Holm equation. Observe that the previous system can be viewed 
as an Hamiltonian system with Hamiltonian function i^(q, p) = | X^j jPiPje"'^'"^-' !. 

In l^j the autors prove the existence of a global multi-peakon solution when strengths pi are positive 
for all i = 1 ... AT; and the convergence of the sequence of multi-peakon solution. If uq is an initial 
data such that the distribution uq — uq^x is a positive Radon measure, there exists a sequence of multi- 
peakons that converges in L°°(M, /^^"^^^(M)). In this case the crucial fact is that no interaction between 
the peakons occurs, and then the gradient remains bounded. However, a general initial data contains 
both positive and negative peakons, as in the so called peakon-antipeakon interaction: one positive 
peakon with strength p, centered in —q, moves forward and one negative anti-peakon in q, with strength 
p moves backward. The evolution of the system produces the overlapping of the two peakons at finite 
time t = T, so that g — > 0. The conservation of the energy E = //(q(t), p(t)) yields 



E = lim — e 



-2kh 



(6) 



and then the quantity p diverges in finite time. At the point (r, 0) occurs thus a singularity for the 
solution u. To extend the solution also after the interaction time with a solution which conserves the 
energy E we can think that at the interaction point emerge an antipeakon/peakon couple, the first, 
negative, moving backward and the second, positive, moving forward with coefficients {—q, —p) and 
{q,p). According to the conservation of the energy, the choice of q and p must satisfy © as t — > r+. It 
yields a change of variables which resolves the singularity at (r, 0) 



c 



■ p^q 



arctan(p) 



with this choice, the Hamiltonian system leads to the ODE 



E 
2 



with 



/(Co;) 



1 



sm' 



and / is a Lipschitz vector field in a neighborhood of the point (■^, |). The solution {({t),uj{t)) of this 
problem provides then the unique couple {q{t),p{t)) which coincides with the classical solution of the 
Hamiltonian system for t < t and extends it for t > t. 

This example suggests the way to construct the multi-peakon solution whenever an interaction 
between peakons occurs (see also .2.^ for an "energetic" motivation). Suppose that two or more peakons 
with strengths pi, ■ ■ ■ ,Pk annihilate at the position q at time r and produce a blow up of the gradient 
Ux- The conservation of the energy yields that there exists and is positive the limit 



lim 



u^{t, x) dx 



where ^~ and are the smahest and the largest characteristic curve passing through the point (t, g). 
Assume that after the interaction appear two peakons with strengths pi, p2 and placed at the position 
gi, 52 • Let consider the change of variables 

z=P2+Pi w = 2arctan(p2 -Pi) r] = q2 + qi C = {p2-Pif{q2-qi), 
then the system © turns out to be 



w 



sm(w) cosh 



C 



2tan2(w/2) 



+ 2z sinh 



C 



2tan2(«)/2) 



j>k+l 



cos\w/2) - sm.\w/2)]e 



isinHsinh^^Ti;^^^^ 



77 

C 



^[l_^g tan2(u,/2)] + 2 COSh 
C 



+ z cosh 



C 

2tan2(u;/2) 



C 



2tan2(to/2) 



j>fc+i 



tani 



j>fc+i 

c 



(«,/2) 



+2C 



sinh 



2tan2(«)/2) 



taii^(w/2) 



+ 



C tan(w)/2) 

- cosh 



+ 



2tan2(«)/2) 



Pi 



Pie 



z cosh 



C 



2tan^(w)/2) 

+ E Pil'i sign {qi - gj)e"l^'"'?^ l 
i>fc+i 



+ tan(t<;/2) sinh 



C 



2tan2(io/2) 



+ 



z cosh 



C 



2tan^(w)/2) 



+ tan(ti;/2) sinh 



C 



2tan2(w)/2) 



+ 



j>k+l 



which is a system of ODE with locally Lipschitz continuous right hand side that can be extended 
smoothly also at the value w = it. The initial data becomes 



(r) = lim y^Pi{t) w{t) = vr rjir) = 2q ({t) = er 



Piir) = lim pi{t) qiir) = lim qi{t) 



k + 1, 



.N 



Thus there exists a unique solution of such a system which provides a multi-peakon solution defined on 
some interval [r, r'[, up to the next interaction time. As in Corollary of [21 Section 7], once we prove 
that Camassa Holm equation is time reversible and the uniqueness of solutions of a Cauchy problem, 
we have that maximal number of peakon interaction is actually k = 2, one with positive strength the 
other with negative one. 



3 A priori bounds 

This section is devoted to the study of some useful properties of the functions u G X^- We start recalling 
an estimate for the L°°— norm of the H^{^) functions. We have 

This estimate give us a bound on the L°°— norm of the conservative solution u of @, in fact the 
conservation of the energy yields 

Mt)\\L^ < \\u{t)\\Hi(R) = Ve^ for every t > 0. (8) 



A 



691: Another fact is the behaviour of the functions u G Xa at the infinitive. If we indicate with C"'^ 
the constant J^{u^ + u^)e"l^l dx, it holds 



supu2(a;)e"l^l <2C°'". (9) 



Indeed, the function 
belongs to //""^(M), moreover 



f{x) = M(t, x)e2 1""'. 



fx = u^e 2 ' ' + ^ sign {x)ue 2 1^1 



and then, by using 0, we have 

1/(^)1' < Wffmm < I + (1 + a)u\'^\y\ dy < . 

Now we study the behaviour at infinitive of the multi-peakon solutions of the Camassa Holm equation. 

Lemma 3.1. (A-priori bounds) Let u be a multi-peakon solution to with initial data u that satisfies 
Then for every t £ M. there exist a continuous function C{t), which depends on C"'" and on the 
energy , such that 

[ [u^it, x) + ul{t, x)]e"l^'l dx < C{t) , (10) 
sup|P,-(i,x)|e"N<C(t). (11) 

||Ma.||Li(R) < C{t) (12) 

Proof. Since |P"| = P", it is sufficient to prove the second inequality with P". Setting 

/(t)= [ [u^{t,x)+ul{t,x)]e'^\''\dx, 



we want to achieve a differential inequality of the kind 

f^m <A+B- m , 

for some constants A and B which depends on the initial data u. We start the discussion proving a 
preliminary estimate for the function P". By definition 



[ P"(t,x)e"l^ldx = J / e"l^'l(ix / e~ 



\^~y\ 



u^{t,y) + 



ui{t,y) 



dy (13) 



from this identity we can use Fubini's theorem to switch the order of the two integrals. Hence we 
compute the following integral 

(' ^a\x\^-\x-y\ ^ _J^f.-\y\ + _J_f.<^\y\ for every ?/ G M . (14) 
Jr 1 - 1 - 

For future use, we observe that the equality ((11]) holds for a E (—1, 1). Substituting (fTl)) in (|13() and 
using the definition of the energy E"^ we have 

/ P"(t, 3;)e"l"ldx < ^. + I{t) . 

Jr (1 -<^V 1 - a 



Having in mind the previous inequality, we are able to estimate the time derivative of the function /. 
From the equations ^ and Q we have 



lit 



I{t) = / [2uut + (u^) J e"l^l dx 



< -2 / (uP" + u<.)x e^l^'l dx < -2u(P" + ui) e' 



+ 2a I |n|(P" + n2)e"l^ldx 



a' 



the previous inequality gives then a bound on the function /, that is 

I{t) < (C7°'^ + ^^/2) exp ( t 



l-a2 



691: To achieve the estimate ((TT|) . set 



i^(t)= P"(t,.)e° 

Proceeding as before, fixed x £ M we compute the derivative w.r.t the time t of the function e~l^' * u^. 



d 



~ I e-\^~y\ul{t,y)dy 
1 



2 ' ^ 



■H 2 



/ e-l^-^lp"(t,y)dy 



< ||'u||LooP"(t,x) + 



|W||L° 

2 

1 „— qIx 



dy 



< ||tx|UcoP«(t,x) + ^e-"l^'l||n||i^i^(i) 
in the same way, the derivative of e"'^' * is 



9 < 



dt 



u{\P^\ + \uUcc\) dy 



<\\u\\l^ r2P"(t,x) + y e-\''-y\p''{t,y)dy 



< \\u\\loo 2P"(t,x) + 



Multiplying the previous 691:two inequalities with e"'^' we get 



d 
'dt 



Kit) < 3 + 



I - a' 



'&K(t) 



which yields 

To achieve the last inequality, we write 



\ux{y)\ dy 



\ux{y)\ dy + 



\ux{y)\ dy 



{y:\u^{y)\e°'\v\ <1} 



{?;:|«,(y)|e-l!'l>l} 



< I e-"l^l dy+ [ ul{y)e''\y\ dy < - + I{t) 
Jr Jr oi 

where the last estimate is given by 



□ 



4 Approximation of the initial data 

In this section we shall construct an approximation with an initial data with a multi-peakon function. 
Our aim is to approximate it with a sequence which has an exponential decay at infinitive uniformly 
w.r.t £. 

Lemma 4.1. Let f G X^. Then for every e > there exists a multi-peakon function g of the form 

N 



g{x) = 5]p,e-l^-*l 



1=1 

such that 



\\f - 9\\m{m) <^ (15) 
/ [g\x)+gl{x)]e''\-Ux<Co (16) 

JR 

for some constant Cq > 0691: which does not depend on e. 
Proof. Let p{x) G Cq" be a cut-off function such that 

• p{x) > 

• p{x) = 1 for every \x\ < 1, p{x) = for every \x\ > 2 

• lRp{x)dx = 1 

and Pe{x) = ^/9(f ) be a moUifiers sequence. Observe that for every e > 0, f{x) = Pe * f{x) is a smooth 
function which approximates the function / in iJ^— norm 

\\f-f\\miR)<Ce (17) 

moreover it belongs to Xa, indeed 



f iPi^) + /'(^)]e"l"l dx < / [ / 691 :m (/^(y) + f^{y))p,ix - y) dy 
JR JR U 

< I [fiy) + fAy)] I Pe{x - y)e"l-l dxdy 

JR JR 

< [ ifiv) + /J(y)]Coe"l^l dy = CoC-'f < oo 

JR 

and Cq is a constant which does not depend on e. From the previous inequality wc can assert that 
for every R > one has ||/||iji(M\[-i?,ij]) < CoCqC""^ uniformly in £ > 0. We can choose thus big 
enough in order to have 

\\f\\miR\l-R.,R.]) < ^/2. (18) 

In the space iJ^([— iig]) we can approximate with a multi-peakon function. By using the 
identity 



the function / can be rewritten in convolution form 
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In the interval [—Re, Re] the previous integral can now be approximated with a Riemann sum 



N 



= ^ Pie 



\x-qi\ 



i=-N 



= ^Re 



f{y)-fxx{y) 



dy . 



Choosing sufficiently large we obtain ||/ — g\\H'^{[~R^,Ri:]) < Together with ((T7)) and ifTH]) this last 
estimate yields the result. □ 



5 Definition of the distance 

In this section we define a metric in order to control the distance between two solution of the equation 
((21) • It is constructed as a Kantorovich-Wasserstein distance. Let T = [0, 27r] the unit circle with the 
end points and vr identified. Consider the metric space (R^ x T,d^), with distance 

d^{{x, u, uj), {x' , u' , uj')) = min{|x — x'\ + \u — x'\ + \uj — uj'\^,,l} 

and for every function u £ Xa, let define the Radon measure on x T 

(7"(^) = / [1 + ul{x)] dx for every Borel set ^ of x T 

J {x(liM.:{x,u{x),2 arctan Ux{x))£A} 

The set of transportation plans consists of the functions ip with the following properties: 

1. ip is absolutely continuous, is increasing and has an absolutely continuous inverse; 

2. sup Ix - V'(a:;)|e°/^'^' < oo; 



3. J^\l-^'ix)\dx< 



oo. 




X 



Figure 1: Transportation plan. 

The conditions [21 and |21 are not restrictive. Indeed, thanks to the exponential decay of functions u,v € 
Xa, the measures a" and a"" located on the graph of u and v respectively, have small mass at the 
infinitive, and then a transportation plan which transports mass from one to the other can be almost the 
identity ip{x) ~ x (see fig. In order to define a distance in the space Xa, we consider an optimization 
problem over all possible transportation plans. Given two functions n, v in X^, we introduce two further 
measurable functions, related to a transportation plan ip: 

(Ai(x) = sup{6'G [0,1] s.t. e-{l + ul{x)) < {l + vl{iP{x)))i;'{x)}, (19) 
4'2ix) = sup {ee[0,l] s.t. l + ul{x) <e- {l + vl{'^{x)))%b'{x)}. (20) 

The functions 4'i,4>2 can be seen as weights that take into account the difference of the masses of the 
measure cr" and a'". In fact, from the definitions (|19j) - (|2n)) one has 

Mx)il + ulix)) = (/.2(V'(x))(l + vl{i;ix)))ij'ix) for a.e. x G M. 
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According to the definitions, the identity max{(/)i(x), (f)2{x)} = 1 holds. Altough the two measures (pia^ 
and have not finite mass, they satisfy (pia^{A) = (l)2(7^{A) for every bounded Borel set ^ C x T. 
Thus, the functions (pi and (p2 represent the percentage of mass actuaUy transported from one measure 
to the other. A distance between the two functions u, v in can be characterized in the following 
way. 

For every ip ^J-^ let X" = {x^u{x),2siTctsia.Ux{x)) and X'' = (^/^(a;), f (-(/'(x)), 2 arctan ^^^(^/'(a;))) and 
consider the functional 

J^{u,v)= [ d^(X",X^)(/>i(x)(l+u2(a;))d2;+ [ \l + ul{x) - {I + vliijix)))^;' {x)\ dx . 
Jr Jr 

Since the previous function is well defined for every ip ^ J-^ we, can define 

Jiu,v) = inf J'^U^v). 

The function J here defined is thus a metric on the space (see [2]). 



6 Comparison with other topologies 

Lemma 6.1. For every u, w E X^ one has 

^ • Ik - AlHR) < J{u,v) < C • ||n - v\\Hi(^^y (21) 

Let (un) be a Cauchy sequence for the distance J such that C"'"" < Cq for every n G N. Then 

i There exists a limit function u £ X^ such that Un ^ u in and the sequence of derivatives Unx 
converges to Ux in IJ'{M) for p G [1,2[. 

ii Let Un the absolutely continuous measure having density w.r.t. Lebesgue measure, then one 
has the weak convergence ^ for some measure fi whose absolutely continuous part has density 

Proof. The first inequality of H21|) can be achieved by estimating the area between the two functions u 
and V. For every ijj £ we can write 



\u — v\dx = I + 
Jr \Jsi Js2 



\u — v\ dx 



where the two subsets Si and ^2 are 



• Si = {x : \x — ij{x)\ < 1 } = Uj[x2j-i,X2j], where in this union we have to take into account that 
these intervals may be either finite or infinite, possibly having xj = ±00 for some j, 

• S2 = {x : \x — Tp{x)\ > 1}. 

The integral over 5*2 can be estimate in the following way: 



\u{x) -v{x)\dx < {\\u\\l-- + \\v\\l--) \x -i;{x)\dx < {E"" + E'')J{u,v). (22) 

S2 Jr 

The last estimate is given by the definition of the functional J. 

As far as the integral over Si is concerned, the integral over Si can be viewed as a sum of the area 
of the regions Aj in the plane M^, bounded by the graph of the curves u, v and by the segments with 



Q 



Figure 2: L — distance between two functions. 



slope ±1 that join the points Quix2j-i) = {x2j~i,u{x2j~i)) and Qv{x2j) = ifp{x2j) , v{ip{x2j))) , where 
{xi} = dSi. We have 

/ |„W-.WM.<y:meas(A,)^ 
is, 

The measure of the subset Aj is the area sweeped by the segment Qu{x) Qv{x). Recahing that in every 
set Aj the function satisfies |x| — 1 < < + 1, a bound on this area is given by 

rx2j 

meas(^j)< / {\x - 'il^{x)\ + \u{x) - v{'il^)\)[l + ul + {1 + vl{ij))ij'] dx 

and then 

/ \u(x) — v{x)\ dx < {\x-i;{x)\ + \u{x)-v{'il))\)[l + ul + {l + vl{ij))ij']dx 
Jsi Jsi 

< J^{u,v) + J^~\u,v) 

this inequahty, together with yields to 

A < C{u,v)J{u,v). 

The proof of the second part of the lemma is perfectly similar to the one of the periodic case, once we 
take into account the exponential decay of the sequence Un- □ 



7 Stability of solutions w.r.t initial data 

Let uq and vq be two multi-peakon initial data. The technique developed in Section El ensures the 
existence of two multi-peakon solution u{t),v{t) for which conserve the energy unless interaction of 
peakon occurs. Suppose then that within a given interval [0, T] no interaction occurs neither for u{t) 
nor for v{t). The aim of this section is to prove the continuity of the functional J w.r.t. the initial data, 
more clearly we prove that there exists a continuous, positive function C{t) such that for t G [0,T] one 
has 

J{u{t),v{t)) < C{t)J{uo,vo). 

Lemma 7.1. If u{t) and v{t) are two multi-peakon solutions defined in the interval [0, T] in which 
no interaction occurs, then there exists a positive, continuous function c{t) which depends only on the 
energies E'^, E"" of the two solutions, such that 

^J{u{t),v{t)) < c{t)J{u{t),v{t)) for all t G [0,r]. (23) 



in 



Proof. We compute the time derivative of the function J'^{u{t),v{t)) with a particular choice of the 
transportation plan ^p = tpf^i^ . Given any ■0o G J^, at every time t G [0, T] we construct '(/'(<) by 
transporting the function ipQ along the characteristic curves. More precisely, since no interaction between 
peakon occurs in the interval [0, T], the functions u(t, •), v(t, •) are Lipschitz continuous, then the flows 
</?^, 99* solutions of the Cauchy problems 

^ifiix) =uit,^iix)) (/p°(x) = x, 

which are the characteristics curves associated to the equation (jlj), are well defined. Now, let x G M. 
'(/'(t) is defined as the composition 

V^(,)(x) = (^t°V'o°K)"'(^)' (24) 

that is 

The function ijj^f-^ belongs to and hence J'^w is well defined, in fact 

1. By the property Q of the function ijjQ and uniqueness of solution of ODE, the function il)^^-^ is an 
increasing function. 

2. Let X E M and ^u{y) be the characteristic curve passing through x at time t. Evaluating \x — 
■0(4) (x)|e"''^l^l along this characteristic curve, and computing the derivative w.r.t. t we obtain 



jy^iy) - ^* (V^o(y))|e"/2|^"(^)l < _ v;(,)(x))| + ||n(t,x)| • |x - ^(,)(x)| 



^a/2\x\ 



by properties ©, and since u, v are Lipschitz continuous in [0,T], there exists two L°° 

functions Ci(t), C2{t) such that 

||x - V'(,)(x)|e-/2|-I < ci(t)|x - Vw(x)|e°/2|-I + c^it) 

by Gronwall Lemma and the hypothesis |x — 'ipQ{x)\e"^'^^^^ < Co, the previous inequality gives the 
property 13 for ipf^^^ 

|x-^(,)(x)|e-/2|-l < Ci{t) = (co + ^*c2(s)ds) elo-^(^)^^ (25) 
3. The last property can be proved by changing the integration variable x = ^uiy) 



|i-^(t)(x)|dx = / \i-i;(^t){xWJ{y)dy= / liviYiy) - {vinMyMiy)\dy 

< / \{^ij{y)-l\dy+ / \{^ly{y)-l\dy+ / \l-i;',{y)\dy. 



Since ^ ^ 

\{^i)'{y)-l\< f \u,{s,x)\-\{ifl)'{y)-l\ds+ f \u,{s,x)\ds 
Jo Jo 

(and a similar estimate for (/?*) and Ux,Vx G by the Gronwall lemma the first two integrals of 
the previous formula are bounded by an absolutely continuous function C{t) in the interval [0, T] 
and then also property 01 holds. 
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At the transportation plan ^(^^ we associate the functions \ (p'^'^ defined according to H19() . (|^n|) . the 
functional J'^w is thus 



dx. 



By deriving J^^^^{u{t),v{t)) w.r.t. t and computing the change of variables along the characteristics, 
the previous derivative can be estimate by the sum of the following terms (we leave out the dependence 
on the integrable variable when it is not essential) 

• h= i Ht, x) - v{t, {x))\4>f (x)(l + ul{t, x)) dx<{l + \\u{t)\\L^ + IIloc ) J^w {u{t),v{t)) , 
P^it, x) - P^{t, (x)(l + ulit, x)) dx , 

'?{l + ulit))dx, 



L 



l + ul{t) 1 + Vl{t,iji^t)) 

the term due to the variation of the base measure 

h = 2 [ dO{lC{t),^''{t))-u^{t){u\t)-P^{t))dx, 



dx . 



and the terms due to the variation of the excess mass 

d 



l + ul{t)-{l+vl{t,i;^t)))i;[^) 



Let us start to estimate the term 12- By definition, the difference of P^ and P'" is written in convolution 
form 

ul{t,yy 



sign {x - y) 



dy 



v\t,i>^t^{y)) + 



_e-|^wW-^(.)(y)l sign(V^(,)(x) - V(t)(y)) 

then in I2 appear the following integrals 

/ {l + ul{t,x)) I e-\''-y\\u\t,y)-v\t,y)\dydx 
Jr Jr 



dx 



A= I u + uZ(t,x]] I e 

B= [{l + ul{t,x)) [ e'\''~y\ sign ix-y)[v\t,y)-v\t,tl^^t)iy))^p'{y)]dy 
Jr Jr 

C= [ {l + ul{t,x)) [ e-l^^-J^I sign (:r-y)-e-l'^w(^)-'^(*)(^)l sign 
Jr Jr 



D = \ I {l + ul{t,x)) 
^ Jr 



v\t,'ip(^t){y)) + -^w-^^'^ xp'^^^iy) dy dx 
l^-^l sign (x - y) [ulit, y) - vl [t, V^(t) (y)] dy 



dx 



A. Switching the order of the two integrals, the term A is bounded by the L^-norm of the difference 
between u and v: 

A < i\\u\\L^ + \\v\\l^) / \u{t,y) -v{t,y)\ / e'l^'^f I (1 + (t, x)) 
< (2 + ^")(||n||Loo + \\v\\L^)\\u{t) - v{t)\\Li 
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and then, by Lemma IHTTl A < C{u,v)J{u{t),v{t)). 
B. Define 



F{y)= r {v\z)-v\^^^){z)Wa){z))dz= [' v\z)dz 



we have, integrating by parts 

/ e~\'''y\s\gTv{x -y)F'{y)dy 
JM. 



<2\Fix)\+ / e-\--y\\Fiy)\dy 
JR 

< \\v\\la.\x - ij{x)\ + / e-\''-y\\y-'ip{y)\dy 



moreover, substituting the previous expression into the term B we obtain 



B< 2EV'^(*)(n(t),t;(t)) + / il+ulit,x)) 



\y - V'(i)(y)l dydx 



2E'^{l^W{n{t),v{t))+ I / {I + ul{t,x))e-\^-y\ dx dy] 



< 2S^(3 + £;") • J^'-'Hu{t),v{t)). 

C. Observe that since the function y i-^ '^{t){y) is non decreasing, the quantities x—y and 

have the same sign, and since the function t ^ e"'*' is Lipschitz continuous either in (— cxd,0) or in 

(0, +00) we have 



-k-s/l _ p-\''P(t)(x)~ip(t)(y)\ 



now 



-min{|x-7/|,|^(t)(x) -^(t)(y)|} < -\x - y\ +2Ci(t), 
where Ci(t) is the function (j^^ . related to the property [21 of -;/'(*) , then 



C< e^iW/ |y-^(y)| 



(i,^{t)(2/)) + 



V'lijly)- / {l + ul{t,x))e-\^-y\dxdy+ 



< 



+2E'' / {l + ui.{t,x))\x -ip(^t){x)\dx 



D. Here we can use the estimate given by the change in base measure. Since 

1 + ul{t,x) - (1 + t;2(t,V'(j)(x)))V(j)(x) dx < J^W(n(t),i;(t)) , 



we obtain 



D < 



I [ {l + nl{t,x)) [ e-l^-^l il + ul{t,y))-{l + vl{t,^^,^{ym[,^{y) 



dy dx 



1 



e J'lsign(x-y) -1 
< {l + E^)J^i^){u{t),v{t))+ [{l + ul{t,x)) 



dy 



< 2(2 + ^«)J'^W(M(t),w(t)) 



dx 



{i'{t){.x)-x) — / {'4'(t){y) -y)dy 



dx 



where in the last estimate we integrated by part as in the term B. 

The control for the terms I3, and /s can be obtained exactly as the ones in [2j, whom we refer 
the reader to. The previous estimates implies that there exists a smooth function C = C^'^{t) which 
depends only to the variable t and to the initial data u, v such that 

j^J^{u,v)<C^'^t)J^{u,v) 

which yields 

J{u{t),v{t)) < J(n(s),T;(s))eli'»*^''''(")H for every s,t£R. 

□ 
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